Fluctuation-dissipation ratio of a spin glass in the aging regime 
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We present the first experimental determination of the time autocorrelation C(t',t) of magneti- 
zation in the non- stationary regime of a spin glass. Quantitative comparison with the response, the 
magnetic susceptibility x(*'> t), is made using a new experimental setup allowing both measurements 
in the same conditions. Clearly, we observe a non-linear fluctuation-dissipation relation between C 
and \, depending weakly on the waiting time t' . Following theoretical developments on mean-field 
models, and lately on short range models, it is predicted that in the limit of long times, the x(C) 
relationship should become independent on t'. A scaling procedure allows us to extrapolate to the 
limit of long waiting times. 

PACS numbers: 



Almost half a century ago, derivation of the fluctu- 
ation dissipation theorem (FDT) jj], ||] which links the 
response function of a system to its time autocorrelation 
function, made it possible to work out dynamics from the 
knowledge of statistical properties at equilibrium. Nev- 
ertheless, this progress was limited by severe restrictions. 
FDT applies only to ergodic systems at equilibrium. Yet, 
such systems represent a very limited part of natural ob- 
jects, and there is now a growing interest on non-ergodic 
systems and on the related challenging problem of the 
existence of fluctuation dissipation (FD) relations valid 
in off-equilibrium situations. 

A way to extend equilibrium concepts to non- 
equilibrium situations is to consider systems in which 
single time dependent quantities (like the average en- 
ergy) are near equilibrium values though quantities which 
depend on two times (like the response to a field) are 
not. Spin glasses 0] are such systems. They remain 
strongly non-stationary even when their rate of energy 
decrease has reached undetectable values. In the ab- 
sence of any external driving force, they slowly evolve 
towards equilibrium, but never reach it, even on geologi- 
cal time-scales. In these conditions, FDT is not expected 
to hold. A quite general FD relation can be written as 
§, D R(t',t) = /3X(t',t)dC(t',t)/dt\ where R(t',t) is 
the impulse response of an observable to its conjugate 
field, C(t',t) the autocorrelation function of the observ- 
able and (3 = l/fe^T. FDT corresponds to X = 1. Deter- 
mination of X, the fluctuation-dissipation ratio (FDR), 
or of an "effective temperature", T e ff = T/X, is the aim 
of many recent theoretical studies which predicted a gen- 
eralization of FDT ||, U in "weak ergodicity breaking" 
systems Q. In the asymptotic limit of large times, it is 
conjectured that the FDR should depend on time only 
through the correlation function: X(t',t) = X(C(t',t)) 
for t' (and t > t') — > oo. The dependence of X on 
C would reflect the level of thermalization of differ- 
ent degrees of freedom within different timescales 
Thus, the integrated forms of the FD relation would 
become x(*'»*) — P $c(t'% X{C)dC (susceptibility func- 



tion) and cr(t',t) = f3j^ {t '' t) X(C)dC (relaxat ion func- 
tion). They would depend on t and t' only through the 
value of C. The field cooled magnetization would read 
Xfc = P io C(M) X(C)dC = /3(1 - ft C(X)dX) (in the 
simplest Ising case with C(t,t) — 1), formally equiva- 
lent to the Gibbs equilibrium susceptibility in the Parisi 
replica symmetry breaking solution for the Sherrington- 
Kirkpatrick model || , with C ^ q (overlap between pure 
states) and X x (repartition of overlap) . Theoretical 
attempts, analytical |9| (with the constraint of stochas- 
tic stability) and numerical jl0| (with the problems of 
size effects), were made in order to confirm the above 
properties in short range models. Up to now, experimen- 
tal investigations correspond only to the quasi-stationary 
regime jll| or are very indirect 

Here we report the result of an investigation of FD re- 
lation in the insulating spin glass CdCr1.7Ino.3S4 |ll|, an 
already very well known compound, with T g = 16.2K. 
Above T g , the susceptibility follows a Curie- Weiss law 
X = C/(T — Q) where C corresponds to ferromagnetic 
clusters of about 50 spins, and 9 w — 9K E^j. The sam- 
ple is a powder with grain sizes around 10 fim, embedded 
in silicon grease to insure good thermal contact between 
grains, and compacted into a coil foil cylindrical sam- 
ple holder 5 mm wide and 40 mm long. The two times 
dependence of the magnetic relaxation (TRM) of this 
compound was extensively studied fl5| . 

In principle, SQUID measurement of magnetic fluc- 
tuations is very simple jl6| 0. The difficulty lies in 
the extreme weakness of the thermodynamic fluctuations 
(of the order of the response to a field about 10 _7 G 
in our case). Therefore, the setup is carefully screened 
against stray fields by superconducting shields, strict 
precautions are taken to suppress spurious drifts of the 
SQUID electronics, and the pick-up coil is a third or- 
der gradiometer. The result is that the proper noise 
power spectrum of the system without sample allows 
time analysis of the magnetic fluctuations signal over up 
to 2000 s of sample fluctuations with more than 20 dB 
of signal/noise ratio. Moreover, in the non-stationary 
regime, the time autocorrelation of magnetic fluctuations 
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C(t',t) = jjYli (m i (t')m i (t)) , where mi is the elemen- 
tary moment at site i, must be determined as an ensemble 
average over a large number of records of the fluctua- 
tion signal, each one initiated by a quench from above 
T g ("birth" of the system). And finally, we want to com- 
pare quantitatively correlation and relaxation data. The 
relaxation function a(t',t) — ( m i(t)) /Hi is mea- 

sured by cooling the sample at time zero from above T g 
to the working temperature in a small field, turning off 
the field at time t and recording the magnetization at 
further times t. Using a classical magnetometer with ho- 
mogeneous field, quantitative comparison between C and 
a is almost impossible due to the strong discrepancy be- 
tween the coupling factors in both experiments. There- 
fore, we have developed a new bridge setup depicted in 
Fig.[l]a, allowing measurements of both fluctuations and 
response. The pick up (PU) coil of self inductance Lq is 
connected to the input coil of a SQUID, of self inductance 
Ls- The whole circuit is superconducting. Relaxation 
measurements use a small coil I inserted in the pick-up 
circuit, and coupled inductively with mutual inductance 
M to an excitation winding. A current Iq injected in 
the excitation results in a field induced by the PU coil 
itself lmG here, clearly in the linear regime though 
inhomogeneous), and the sample response is measured 
by the SQUID. To get rid of the term L , the sample 
branch is balanced by a similar one without sample, ex- 
cited oppositely (see Fig.|l|a). The flux delivered to the 
PU by an elementary moment m at position r is given 
by <&(i) = m(r, i)h(r) where h is the magnetic field pro- 
duced by a unit of current flowing in the PU. Flux con- 
servation in the PU circuit results in a current Is flowing 
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FIG. 1: a) Schematic of the detection circuit. The pick-up 
coil (right side), containing the cylindrical sample, is a third 
order gradiometer made of +3 -6 +6 -3 turns, b) Calibration 
is obtained by measuring relaxation versus correlation in a 
high conductivity copper sample at equilibrium at 4.2K . 



in the input coil of the SQUID whose output voltage is 
Vs = GIs- Detailed analysis of the system will be pub- 
lished elsewere. The main features are as follows. 

As the fluctuations of elementary moments in the sam- 
ple are homogeneous and spatially uncorrelated at the 
scale of the PU, the SQUID output voltage autocorrela- 
tion is given by: 



Cs(t',t) = {V s (t')Vs(t)) = C(t',t)Q 
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Q = h 2 (ri) where the index i refers to a moment site, 
is the coupling factor to the PU, including demagnetizing 
field effects since h is the internal field. 

The elementary moment response at site i is Ri(t', t) 
a (mi(t)) /9h(r,-, t'). Taking into account that the medium 
is homogeneous, the relaxation function of the SQUID 
output voltage is given by 



a s (t',t) = ^^- = cr(t',t)Q 
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Thus, the coupling factor Q disappears in the relation be- 
tween Cs and as, independently on the nature and shape 
of the sample. There remains only the inductance terms 
M, Lq and Ls- These being difficult to determine with 
enough accuracy, absolute calibration was performed us- 
ing a copper sample of high conductivity, by measur- 
ing as(t',t) and Cs(t',t) — computed by standard fast 
Fourier transform algorithm — at 4.2K ( 4 Hc boiling tem- 
perature at normal pressure): with this ergodic material, 
the relation between both measured quantities is linear 
with slope A/T, where A is the sample independent 
calibration factor (see Figj^b). From the knowledge of A, 
determined at 4.2K, the system is equivalent to a ther- 
mometer, i.e. the FDT slope is known exactly at any 
temperature. 

In the spin glass sample, Cs(t',t) and as(t',t) were 
measured at T = 0.8T g after quench from a temperature 
T « 1.2Tj. To get a precise definition of the "birth" 
time, a minimum value of 100 s was chosen for t' . The 
autocorrelation was determined from an ensemble of 320 
records of 12000 s of the fluctuation signal. The ensemble 
averages were computed in each record from the signal at 
t', averaged over St' < t'/20, and the one at t, averaged 
over St ^ (i — t')/10 — the best compromise allowing a 
good average convergence still being compatible with the 
non-stationarity — , and averaging over all records. As 
there is an arbitrary offset in the SQUID signal, the con- 
nected correlation was computed. Nevertheless, this was 
not enough to suppress the effect of spurious fluctuation 
modes of period much longer than 2000 s, giving a non- 
zero average offset on the correlation results. Thus, as a 
first step, we have plotted all correlation data, taking as 
the origin the value of (Vg(t')). Due to the elementary 
measurement time constant this last term corresponds 
to an average over t — t' about I0~ 2 s, i.e. a range of 
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(t — t')/t' corresponding to stationary regime. Thus, all 
Cs data are shifted by a common offset Co . The result is 
shown in Fig.|^a (right sided scale), as a function of t — t' 
for values of t' from 100 s to 10000 s. Residual oscillations 
— and large error bars — for t' — 100s reveal the limit of 
efficiency of our averaging procedure. Corresponding re- 
laxation data are plotted on FigJ^b. In both results, one 
can see that the curves merge at low t — t' , meaning 
that they do not depend on t' (stationary regime). At 
t — t' ^ t', they strongly depend on t' , the slower decay 
corresponding to the longer t' . 

The correlation offset must be determined. As zero 
of correlation is unreachable in experimental time, cor- 
rection of the offset could be obtained from the knowl- 
edge of C(t,t). Nevertheless, due to clustering, C(t,t) 
depends on temperature and cannot be determined from 
the high temperature susceptibility. In canonical com- 
pounds like 1% Cu:Mn ||, with negligible clustering, the 
field-cooled susceptibility is temperature independent in 
agreement with the Parisi- Toulouse hypothesis jl0|, p^[ , 
yielding C(t,t) = T q xfc{T)- We used a generalization 
of this relation with the condition that a smooth depen- 
dence of C(t, t; T)/T must result [[Hi. This was obtained 



by using for T g a slightly different value, T* 



Then, from the value of the calibration factor 
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FIG. 2: Aging and scaling of (a) correlation (b) relaxation 
at T — 0.8T 9 . Both are measured for waiting times t' = 100 
(o), 200 (A) , 500(V), 1000(o), 2000 (+) , 5000(x), 10000(*) 
seconds from bottom to top. Reported error-bars on correla- 
tion have a length of two standard-deviation, corresponding 
to averages over records. In insets, scaling of the aging parts 
versus < = (i 1_M - f' 1_M )/(l - fx), using fx = 0.87. The sta- 
tionary parts are found to obey a power-law decrease with an 
exponent a — 0.05. 



17.2K. 

A, and 

writing C(t,t;0.8T g ) = 17.2 XF c(0.ST g ), C s (t, t; 0.8T g ) 
can be determined, and suppression of the offset can 
be performed by using the x(C) plot, first introduced 
by Cugliandolo and Kurchan 0. We plot the normal- 
ized susceptibility function x(i',t) = 1 — cr(t',t) where 
a(t',t) — <rs{t' ,t)/as(t,t) (note that <r(t,t) = xfc) ver- 
sus normalized autocorrelation C(t', t) — Co = (Cs(t',t) — 
Co)/Cs(t,t;T) for all experimental values of t' . In this 
graph, the FDT line has slope —T*/T and crosses the 

C axis at C = 1. On the data, a clear linear range ap- 
pears at large C (small t — t'), displaying the FDT slope 
with error < 3% in the sector C ^ 0.47. This allows the 
suppression of the correlation offset by horizontal shift of 
the data. The result is shown in Fig.|[ It is of course 
based on a rough ansatz on C(t',t;T) which needs fur- 
ther justifications, but we stress that the induced uncer- 
tainty concerns only the position of the zero on the C 
axis, and not the shape and slope of the curves. With 
decreasing C (increasing t — t'^t'), the data points de- 
part from the FDT line. Indeed, C(t',t —> oo) = and 
Xfc = x[t', t — > oo) = 1. The mean slope of the off FDT 
data corresponds to a temperature of about 30K. This 
value is far above our annealing temperature, ruling out 
a simple interpretation in terms of a "fictive" tempera- 
ture pcfl . Despite the scatter of the results, a tendency 
for the data at small t' to depart the FDT line at larger 
values of C is clear: it is experimentally impossible to 
fulfill the condition of timescales separation underlying 
the existing theories . Even if the long t' limit for x(C) 
does exist, it is not reached iii_the plot of data in Fig.|| 
and a t' dependence of the x{C) curves is expected. 
^ The left sided scales in Fig.^Ja and b correspond to 
C(t',t) and a(t',t) respectively. In former works, it was 
shown that the whole relaxation curves could be scaled 
as the sum of two contributions, one stationary and one 
non-stationary jL5| 

a(t',t) = (1 - A)(l + (t- t')/t )- a + Ap(C), (3) 

where to is an elementary time of order 10~ n s, ip is 
a scaling function of an effective time parameter £ cx 
— t' 1 ^^ depending on the sub-aging coefficient fi < 1 
Jlo| , and a can be determined with good precision from 
the stationary power spectrum of fluctuations S(u>) cx 
w q-i mse t m Fig.||b displays the result of the scal- 
ing on the relaxation curves with a = 0.05, A = 0.21 
and fx — 0.87. As shown in the inset of Fig.^ja, the scal- 
ing works rather well on the autocorrelation curves with 
the same exponents, but now, qEA, the Edwards Ander- 
son order parameter, replaces A. We get qsA = 0.37. 
These results show clearly that the stationary part of 
the dynamics is still important yet in the aging regime, 
i.e. that the limit of long t' is not reached within the 
timescale of our experiments (in fact, timescale separa- 
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FIG. 3: FD-plot. Relaxation measurements are plotted versus 
correlation functions for each t'. The dot-dashed line (FDT 
line) is calculated for T = 0.8Tg — 13. 3K, from the cali- 
bration obtained with the copper sample. The dashed line 
represents the scaling extrapolation for t' —* oo. The branch- 
ing point with the FDT line, corresponds to C = Qea (square 
symbol, with size giving the error range). In Inset, the same 
data in the whole range. 



tion is realized if t' ^ r where r is the observation time 
such that C s tat{r) < qsA)- 

If granted, the scaling gives the long time limit of 
the non-stationary part of the dynamics, allowing a plot 
of the long times asymptotic non-stationary part of the 
x(C) curve. Of course, here we verify it only over 2 
decades of time, up to t' = 10000 s, but it was proven to 
be relevant on TRM up to t' = 100000 s @. The dashed 
line in Fig.^ is obtained by plotting the smoothed curves 
of aging parts of x(C) versus C{Q). According to theo- 
retical conjectures, dx(C)/dC would represent the static 
quantity x(q) One can see that the curve does not 
point exactly towards x = 1 but about 5% below. There- 
fore i) either the ansatz used to determine C(t, t; T) is not 
realistic enough, ii) or the time scaling is no longer valid 
at the very large t 1 needed for timescales separation. For 
future progress, scaling developments outside the strict 
time range separations at the basis of the "adiabatic cool- 
ing" analysis or the "weak memory" analysis [^) are 
needed. It seems that such developments are out of the 
theoretical possibilities for the moment. Toy models |||] 
presently under development could allow a phenomeno- 
logical approach of the problem. 

In conclusion, we have presented the first experimental 
determination of the non stationary time autocorrelation 
of magnetization in a spin glass, an archetype of a com- 
plex system. With the help of the time scaling properties 
of both the relaxation and the autocorrelation, we were 



able to propose a first experimental approach of a pos- 
sible generalization of FDT to non-stationary systems. 
Results at several temperatures are now needed in order 
to get a complete description of the x(C) behavior in the 
whole temperature range. 
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